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a modified  Eward  solution  forms  the  basis  for  the  planar  wing  cases  and 
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experimental  data.  The  gapped  elevon  cases  for  supersonic  flow  is  also 
analyzed  and  preliminary  results  have  identified  thickness  effect,  as  being 
very  important,  at  least  for  small  elevon  deflections. 
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SUMMARY 


The  aerodynamic  loading  for  deformed  wings  in  both  subsonic  and 
supersonic  flow  has  been  under  study  for  the  past  year.  The  basic  solu- 
tion technique  falls  into  the  potential  flow  category  with  appropriate 
restrictions.  A lifting  surface  kernel  function  formulation  is  used  for 
both  subsonic  and  supersonic  flow  and  results  are  obtained  which  agree 
very  well  with  experimental  data.  Subsonic  solutions  for  deformed  wings 
with  deflected  elevons  are  obtained  in  a semi-closed  form  summation 
manner.  Cases  under  study  which  include  gaps  between  the  elevon  and  wing 
are  as  yet  incomplete  but  is  presently  being  pursued.  For  supersonic  flow, 
a modified  Eward  solution  forms  the  basis  for  the  planar  wing  cases  and 
vorticity  paneling  is  added  to  account  for  local  deformations  in  the  wing 
mean  camber  surface.  Results  are  obtained  which  agree  very  well  with 
experimental  data.  The  gapped  elevon  cases  for  supersonic  flow  is  also 
analyzed  and  preliminary  results  have  identified  thickness  effect,  as  being 
very  important,  at  least  for  small  elevon  deflections. 
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I . INTRODUCTION 

The  aerodynamic  loading  for  deformed  wings  in  both  subsonic  and 
supersonic  flow  has  been  under  study  for  the  past  year.  The  problem  of 
wing  loading  solutions  have  generally  been  confined  to  potential  flow 
formulations.  These  solutions  have  been  applied  primarily  to  planar 
wings  but  in  some  cases,  such  as  vortex  lattice,  are  applicable  to  the 
deformed  wing  problem.  The  present  formulation  falls  into  the  potential 
flow  category,  but  unlike  vortex  lattice  techniques,  is  of  the  continuous 
vorticity  distribution  type  yielding  continuous  loading  profiles. 

Not  only  are  local  wing  deformations  included  in  the  analysis  but 
for  subsonic  flow,  the  elevon  loading  is  computed  for  displaced  elevons 
(gap  between  elevon  and  wing  trailing  edge).  In  all  cases,  sectional 
as  well  as  total  loads  and  moments  are  computed  in  a "semi"  closed  form 
manner.  A comprehensive  discussion  of  the  basic  formulation  for  subsonic 
flow  is  presented  in  Ref.  1. 

For  supersonic  flow,  a similar  kernel  function  approach  is  taken 
leading  to  general  functional  forms  for  the  localized  loading  on  wings 
in  supersonic  flow.  A more  complete  derivation  of  the  applicable 
equations  for  supersonic  flow  is  presented  in  Ref.  1.  In  orde.  to 
preserve  continuity,  however,  the  general  equations  applicable  for  both 
subsonic  and  supersonic  flow  regimes  are  repeated  in  this  report. 
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II.  SUBSONIC  FLOW 


Fundamental  Equations 

For  a planar  lifting  surface,  the  potential  equation  is 


<t>(x,y,z)  = 


i_  r t , /. 

4’  (y-yc)2+z2  y 


f(x-x  )2+B2(y-y 
o o 


— ) 

-y  )2+B2z2/ 


dx  dy  ( 
o o 


The  closed  form  integration  technique  for  wing  sub-elements  as  described 
in  Ref.  1 yields  excellent  results  for  planar  wings  with  the  assumed 
pressure  loading 


N M rz — T 

ACp(C,n)  = l l BL  sin(2nH-l)ey^  C 

n=0  m=0 


Spanwise  Chordwise 


The  applicable  boundary  condition  is 


w(x,y)  + sin  a(x,y)  = 0 


and  solutions  are  obtained  for  the  various  sectional  coefficients  as 


cc.(y  ) = y E sin(2m+l)9 
x,  o _ m 

m=0 


ccs(y  ) M sin(2m+l)0 

ccd(yo)  = 8 ^ (2m+1)Em  -JUT 1 ^ 

m=0  o 


M N 

cc  (y  ) * -c(y  ) 7 T B I . . sin(2m+l)9 

mwo  o L n+1  o 

m=0  n=0 


where  L = Nm  + n. 


H 

I 


i 
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Total  forces  and  moments  are  then  found  by  integration  as 


IT 

16S 


M 


I 

m=0 


(2nri-l)E  2 
m 


(7) 

(8) 


m 


1 M 

— [<c_-c  ) l F G 
Sc  R T m=0  m m 


M 


' 2 CR  F(1)  'tan  A ^ E™  G»] 


m=0 


m m 


(9) 


Comparison  with  experimental  data  is  excellent  (Ref.  1)  for  a wide  variety 
of  wing  shapes. 

Order  of  Polynomials 

Before  proceeding,  comments  on  the  effect  of  polynomial  order  on  the 

solution  is  in  order.  The  unknown  constants,  B in  Eq.  (2)  are  really 

J-i 

coefficients  for  the  various  terms  in  the  polynomial  describing  the 
pressure  loading  on  the  wing.  It  has  been  suggested  that  if  the  polynomial 
order  is  increased  to  a large  number  of  terms  (i.e.,  > 50)  that  the 
sectional  coefficients  such  as  cc^  will  tend  to  "oscillate"  in  the  span- 
wise  direction.  For  the  formulation  used  in  the  present  analysis,  no 
oscillations  are  observable  at  least  up  to  50  terms  in  the  equation.  How- 
ever, oscillations  do  indeed  occur  if  terms  on  the  order  of  80  or  more  are 
used  in  the  pressure  loading.  It  does  not  appear  that  this  is  due  to 
"numerical"  instability  as  would  be  the  case  if  numerical  Integration  were 
used,  but  rather  due  to  overkill.  That  is,  the  pressure  loading  function 
is  rather  general  and  for  most  wings  only  a few  terms  (<  20)  are  required 
for  a solution.  If  more  terms  are  used,  it  is  analogous  to  fitting  a 50th 
degree  polynomial  to  a linear  curve;  oscillations  are  the  inevitable  results. 
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Elevon  Loading  and  Hinge  Moments 

For  the  past  year,  concentration  has  been  placed  on  predicting  the 
loading  of  a full  span  elevon  placed  in  the  wake  of  a lifting  surface  and 
on  the  associated  elevon  hinge  moments.  Displacement  of  the  elevon  creating 
a wing-elevon  gap  and  its  effect  on  the  wing-elevon  loading  has  also  been 
under  study.  It  is  apparent  that  when  a gap  exists  between  the  wing  and 
the  elevon,  the  loading  distributions  on  each  will  change.  For  a gap,  the 
wing  and  elevon  are  treated  as  two  wings  with  the  requirement  that  the  gap 
is  less  than  local  wing  chord,  and  therefore,  the  wing  wake  has  no  signifi- 
cant rollup.  Two  loading  functions  are  used,  satisfying  the  infinite 
pressure  differential  at  the  leading  edge  and  the  Kutta  condition  at  the 
trailing  edge  of  both  the  wing  and  the  elevon.  The  assumed  loading  for 
both  the  gap  and  the  gap-sealed  cases  are  shown  schematically  in  Figure  1. 

The  sectional  lift  for  a wing  with  an  elevon  is 

x=c  x~C 

cc  (y  ) = / w AC  (5,n)dx  + / To  AC  (£,n)dx  (10) 

° x=0  Pw  x=(c  +e)  PE 

w 

where  the  total  chord  length  is 


c„  = c + c„  + e 
T w E 
o 


(ID 


and  e denotes  the  gap  width. 

For  the  wing,  the  nondimensional  chordwise  variable  is  defined  as 


= (x-Xtt,  >/c. 


w 


LE 


w 


w 


(12) 


J 
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and  the  sectional  lift  for  the  wing  becomes 


cc  - c / AC  <C  »n)  d£ 
A w ' p w 
w C«0  w 


Similarly,  for  the  elevon, 


5e  ■ <*--LE  )/cE 
£ 


where  the  elevon  leading  edge  is  now 


XLEg  " *L^+  Cw  + E (15) 

The  sectional  lift  is  now  given  by 

?w=1  5E=1 

cc  (y  ) = c / AC  U,n)  d£;  + c / AC  (S,n)  dC  0-6) 

AO  W.  -P  - p_ 

Cw-°  Pw  CE=0  PE 

Each  integral  is  treated  separately  and  the  integration  is  carried  out  as 
with  the  wing  alone  analysis  of  Ref.  1.  The  total  lift  and  moment  is  also 
similarly  found.  It  should  be  noted  that  the  gap  dimensions  are  small 
and  are  excluded  in  defining  all  reference  dimensions.  Thus  the  total 


wing  area  is 


S = S + S_ 
w E 


Extension  of  this  analysis  to  include  hinge  moment  computations  was 


also  done  in  Ref.  1.  However,  by  definition  the  hinge  line  was  located 
at  the  leading  edge  of  the  elevon.  For  other  hinge  line  locations,  the 


moment  coefficient  becomes 


c = c + c (x^  /c  ) 
“HL  mLE  LE  ^ E 
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where  is  the  chordwise  location  of  the  hinge  line  defined  by 

x,,T  = x - x,„  -c  -e.  It  has  been  shown  that  the  sectional  lift  for  the 


*LE  "Cw 
w 


elevon  is 


where 


nE  I mE 

cc,  - CB  l ~ l IB  sin  (2q+l)  0 

*E  E p-0  CE  q-0  P L 


I = n/2 
o 


I = (f2^-)  I - 

p 2p+2  p-1 


L = (Ne)  q + p. 


The  total  lift  produced  by  the  elevon  is 


CT  = / cc  (n)dn 

E E 0 E 


Using  (19)  in  Eq.  (23)  yields  (after  integration) 


where 


“L  2S 
E E 


Ne 

H = l Bt  1 

q p=0  L P 


In  a similar  manner  the  moment  produced  about  the  elevon  hinge  line  is 


C = C + 

“HL  mLE  2SE  CE 


where  C and  C are  moments  produced  by  the  elevon  only. 

DhL  \e 
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Wing  Deformations 


A wing  which  is  warped  such  that  all  points  on  the  wing  do  not  lie 
in  the  x-y  plane  presents  a different  view  to  the  flow  field  from  that 
of  the  planar  wing  considered  previously.  To  account  for  wing  warpage, 
it  is  necessary  to  consider  the  angle  that  each  control  point  and  each 
grid  element  makes  with  the  x-y  plane. 

For  warpage  in  the  chordwise  direction,  the  downwash  produced  by 
each  grid  element  is  proportional  to  the  angle  which  that  element  makes 
with  the  x-axis,  or 


Aw(x’y)warped  = [Aw(x’y) unwarped1  cos  “i  < 

where  is  the  wing  slope  angle  relative  to  the  x-axis  and  is  given  by 


= tan-1  (dz/dx)^ 


The  boundary  condition  is  for  no  flow  through  the  wing,  or  that  the  normal 
component  of  velocity  to  the  wing  be  zero.  The  angle  that  the  control 
point  makes  with  the  x-axis  is  given  by 


a = tan-1(dz/dx) 
c c 


so  that  the  downwash  from  a sublement  at  a control  point  is 


Aw(x,y)  . - [Aw(x,y)  .]  cos  (a  -a  ) (29 

warped  unwarped  i c 

Definition  of  a wing  deformation  function  z(x)  of  the  mean  camber  line 
can  therefore  produce  many  variations  and  amplitudes  of  wing  distortions 
in  the  chordwise  direction.  However,  the  condition  of  no  flow  separation 
still  applies  at  all  points. 


» 

M 

» * 

U 

I 


! 

i 


9 


Subsonic  Results 

It  has  been  shown  (Ref.  1)  that  the  preceeding  wing-elevon  analysis 
produces  results  which  agree  well  with  experimental  data.  However,  the 
availability  of  experimental  data  for  a gapped  wing-elevon  combination, 
or  for  deformed  wings  is  very  limited;  therefore,  for  the  present  time, 
the  investigation  is  restricted  to  predicting  results  and  trends  of  the 
theoretical  analysis. 

The  research  was  conducted  using  two  types  of  wings,  a rectangular 
wing,  and  a low  aspect  ratio  delta  wing.  Four  variations  were  used  for  each 
wing  planform;  a planar  chord,  a wing  with  twist,  a wing  with  circular  arc 
camber,  and  a wing  with  a sinusoidally  deformed  camber  line.  Hinge  moments, 
total  moment,  and  total  lift  were  then  computed  for  varying  gap  sizes  and 
angles-of-attack . 

The  computations  for  the  gapped  elevon  cases  are  as  yet  incomplete 
and  will  therefore  not  be  presented  in  this  report.  At  present,  it  may  be 
stated  that  gaps  at  the  wing-elevon  juncture  certainly  change  the  loading 
on  both  the  wing  and  elevon.  Quantitative  calculations  regarding  this 
effect  will  be  presented  in  a later  report. 

For  the  non-gapped  case.  Fig.  2 illustrates  C and  C for  a planar 

n H f. 

a 6 

versus  a sinusoidally  deformed  wing.  The  hinge  moment  slopes  for  a deformed 
wing  is  essentially  the  same  as  the  planar  wing,  differing  only  by  a 
constant.  It  is  seen  also  that  the  hinge  moment  decreased  linearly  with 
increased  leading  edge  camber  angle. 

Several  other  configurations  are  currently  being  investigated  with 


varying  types  of  deformations  as  listed  in  the  following  table. 


TABLE 


1.  Rectangular,  planar  wing. 

2.  Rectangular  wing,  6°  twist  (washout)  at  the  tip. 

3.  Rectangular  wing,  circular  arc  camber 

c u 

z(x)  » [r2  -(x  — ^)2]* 

4.  Rectangular  wing,  sinusoidal  deformation 
z(x)  ■ A sin  (2ttx) 

5.  Delta,  planar  wing 

6.  Delta  wing,  6°  twist  (washout)  at  the  tip 

7.  Delta  wing,  circular  arc  camber 

z(x)  - [r2  -(x  - 

8.  Delta  wing,  sinusoidal  deformation 
z(x)  * A sin  (2irx) 

For  each  configuration,  the  dependence  of  the  hinge  moment  coefficient 
on  the  hingeline  location  and  on  the  gap  size  is  to  be  determined. 
Dependence  on  Mach  number  and  angle  of  attack  is  still  under  investigation 
with  results  expected  in  the  near  future.  Definitive  efforts  are  also 
under  way  to  quantify  the  dependence  of  gap  size,  etc.,  on  total  lift  and 
pitching  moment  of  the  entire  configuration. 


III.  SUPERSONIC  FLOW 


| 


i 


The  analysis  used  in  the  present  research  effort  for  wings  in 
supersonic  flow  Is  based  in  part  on  Ewards  theory  and  in  part  on 
vorticity  paneling  leading  to  simple  expressions  for  the  basic  aero- 
dynamic coefficients.  As  in  subsonic  flow,  the  wing  is  subdivided  into 
a number  of  small  subpanels  over  which  it  is  assumed  that  the  panel 
strength  is  constant.  The  vorticity  panel  is  actually  a continuously 
distributed  doublet  panel  which  satisfies  the  supersonic  potential 
equation 


dx  dy 
o o 


(30) 


. y(x  ,y  ) / (x-x^ ) \ 

4>(x,y,z)  = y-  jf  2 — 2 z j — - | 

Sq  (y-yQ)2+z2  1 y (x-xq) 2-e2 (y-yD) 2~B2z2 J 

As  for  subsonic  flow,  the  vorticity  distribution  Y(x0»yQ)  n^y  be  replaced 
by  the  pressure  loading  coefficient  ^AC^x^jy^)  • The  downwash  equation 
is  obtained  by  differentiation  in  the  standard  manner  and  when  evaluated 
in  the  z=0  plane  (planar  flow)  becomes 


. AC(x.y)  (x-x) 

w(x  v)  _ A-  f / P-_g_  2_  ° 

4ir  sq  (y-yQ)2  ^( x-Xq ) 2- 0 2 (y-yQ ) 2 


dx  dy 
o o 


(31) 


where  the  area  S is  that  contained  on  the  wing  surface  in  the  upstream 
o 

running  Mach  cone  emanating  from  the  field  point  (x,y). 

The  assumed  pressure  loading  differs  from  that  previously  reported 
in  Ref.  1 and  is  written  as  the  sum  of  two  expressions  as 


N M 


AC  (x  ,y  ) = P(x  ,y  )sin(a+6  ) + P(x  ,y  ) \ J B n”1  cos 


p o o 


o o 


0"0  “ L 

n=*0  m=0 


(32) 


. 


12 


13 


The  first  term  on  the  right-hand  side  is  a modified  Eward  solution 

to  account  for  local  wing  deformations  and  is  valid  for  both  subsonic 

and  supersonic  leading  edges.  The  second  term  uses  the  pressure  functions 

P(x  ,y  )'  as  a weighting  function  which,  when  multiplied  by  the  unxnown 
o o 

coefficient  B , and  the  series  n™  cos  accounts  for  small  deformations 

in  the  wing  surface  in  the  upstream  Mach  cone  from  the  point  (x,y).  The 

order  of  the  terms  in  the  second  expression  seldom  exceeds  3 or  4 and 

for  planar  wings,  the  entire  second  term  is  unnecessary.  Increasing  the 

order  for  larger  wing  deformations  or  higher  angles  of  attack  does  not 

seem  to  provide  better  agreement  with  experimental  data. 

As  in  subsonic  flow,  the  AC  term  is  removed  from  the  integral  in 

P 

Eq.  (31)  and  the  resultant  expression  is  integrated  in  closed  form  over 
a small  subpanel.  Results  of  this  integration  over  various  shaped  sub- 
elements is  presented  in  Ref.  1. 

It  is  convenient  to  define  the  value  of  this  integral  as  3AK  so  that 
the  total  downwash  at  a point  (x,y)  becomes 


w(x,y) 


i=l 


(AK)  [AC  (x  ,y  ) J 
i pool 


(33) 


where  I is  the  total  number  of  subpanels  in  the  region  (cone)  of  integra- 
tion. The  boundary  condition  is 

w(x,y)  + sin  (a+6  ) = 0 (34) 

c 

where  0 is  the  local  deformation  angle  at  a control  point.  Rearranging 
c 

and  combining  Equations  (32),  (33)  and  (34)  yields 
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1 

l (AK)  [P(x  ,y  ) sin  (a+6  )]  + 

, i o o o l 


I 


I 

1=1 


(AK).[P(x  ,y  ) 
i o o 


N M 


l l 

n=0  m=0 


m nirC, 
ri  cos  — y~ 


= ~ f sin(a+6c) 


(35) 


where  as  before  L = n + mN.  The  above  expression,  Eq.  (35),  is  con- 
veniently arranged  in  matrix  form  for  a Gaussian  reduction  solution  of 
the  unknown  coefficient,  B^. 

Body  volume  effects  on  the  wing  loading  are  accounted  for  by 
assuming  that  the  body  is  cylindrically  shaped.  The  cylinder  is  generated 
with  an  infinite  line  doublet  whose  axis  is  in  the  z-direction  perpen- 
dicular to  the  wing  planform.  (Lifting  interference  effects  are  not 
accounted  for  in  the  present  analysis  but  is  under  study  for  inclusion 
in  future  work. ) The  result  of  adding  the  body  volume  to  the  analysis 
simply  adds  another  term  to  the  right-hand  side  of  Eq.  (35)  which  after 
a little  manipulation  becomes 

4ira2  sin  a 
8 (a+n) 2 

where  a is  the  nondimensional  body  radius  defined  by  a = r^/(b/2). 

Unlike  the  subsonic  analysis,  the  (supersonic)  sectional  properties 
cannot  be  integrated  in  closed  form  but  must  be  done  numerically. 

Results  for  Supersonic  Flow  - Planar  Wings 

Several  wing  planform  shapes  were  analyzed  for  which  both  undeformed 
and  deformed  data  were  available.  A sketch  of  the  four  planforms  are 


fl 
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shown  in  Fig.  3.  Results  of  the  pressure  loading  on  these  wings  in  an 
undeformed  state  are  presented  in  Figs.  4-9  for  Mach  numbers  of  1.61  and 
2.01.  In  these  figures,  the  theoretical  solution  is  that  taken  from 
Refs.  2 and  3.  Minor  improvements  could  possibly  be  made  using  vorticity 
paneling;  however,  since  agreement  with  experimental  data  is  excellent, 
this  is  unnecessary.  Experimental  data  for  both  the  planar  and  deformed 
wings  are  obtained  from  Refs.  4 and  5. 

It  is  apparent  from  Figs.  4 through  9 that  the  linear  theory  is 
excellent  for  subsonic  and  supersonic  leading  edges.  Th°  disagreement  near 
the  leading  edge  and  wing  tip  is  to  be  expected  and  is  due  to  wing  thick- 
ness, shocks  and  wing  twist;  the  latter,  of  course,  is  not  significant 
for  the  low  aspect  ratio  wings. 

In  Figs.  8 and  9 , the  error  is  attributed  primarily  to  wing  twist 

where  the  wings  under  wind  tunnel  tests  experienced  up  to  1.5°  aeroelastic 

twist  (washout)  at  the  wing  tips.  Washout  in  the  vacinity  of  the  tip 

would  tend  to  lower  the  AC  distribution  which  is  indicated  in  the  experi- 

P 

mental  results. 

Deformed  Wings 

These  same  wings  under  various  kinds  of  deformation  were  also 
analyzed  by  ineluding  vorticity  paneling  as  well  as  the  previous  theoretical 
loading  terms  in  the  ACr  distribution.  Deformation  shapes  of  three  of 
these  wings  are  shown  in  Fig.  10.  The  delta  wing  (#1)  in  Fig.  3 was 
cambered  with  a mean  camber  line  curve  fitted  and  defined  by 

(z/cR)  = tan  3° [•£■  - (n+£(l-n)-.5)2] 


(36) 
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Fie.  3.  Dimensional  Sketch  of  Four  Wing  Planform  Shapes  Under  Study. 
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g.  4.  Pressure  Distribution  on  a Planar  Trapezoidal 
Wing,  M-1.61. 
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Fig.  10.  Schematic  of  Three  Deformed  Wing  Shapes. 


24 


i 

i 

i 


where  n**y/(b/2)  and  £*(x-x^g)/c.  The  camber  (See  Fig.  10)  is  uniform  in 
the  spanwise  direction  and  produces  a maximum  deflection  angle  of  three 
degrees  relative  to  the  chord  line.  Although  this  seems  like  a small 
task,  this  is  actually  a severe  test  of  the  theory  and  its  agreement  with 
experimental  data.  In  Figs.  11  and  12,  results  are  presented  for  the 
cambered  delta  wing  at  Mach  numbers  of  1.61  and  2.01.  Agreement  at  the 
inboard  stations  are  relatively  good  while  outboard  agreement  is  rather 
poor.  However,  note  that  in  the  regions  where  most  of  the  wing 
loading  occurs,  the  agreement  is  satisfactory.  Additionally,  the  inclusion 
of  the  doublet  paneling  terms  in  the  mathematical  formulation  provides  the 
correct  trends  and  an  added  degree  of  accuracy. 

For  the  same  planform  shape,  linear  spanwise  twist  is  added  to  the 
camber  deformation  such  that  the  mean  camber  surface  is  now  defined  as 


(z/cR)  = - .5 (tan  8°) (n) (n+£(l-n))2  (37) 

Results  of  this  configuration  are  shown  in  Figs.  13  and  14. 

Similar  results  are  obtained  for  these  cases  as  was  observed  for  the 
cambered  (alone)  case.  Agreement  is  acceptable  at  the  inboard  stations 

j 

but  is  questionable  at  the  outboard  stations. 


In  observing  the  experimental  data  for  both  the  cambered  alone  and 
the  cambered  and  twisted  wings,  the  outboard  stations  (n  > -5)  appear  to 
have  a supersonic  leading  edge.  Whereas  the  high  leading  edge  sweep  angle 
gives  rise  to  a subsonic  leading  edge.  In  addition,  at  the  n = 0.5  station, 
a weak  shock  appears  to  be  located  near  C = «4  for  M = 1.61  but  is  not 
present  for  M = 2.01.  Since  the  theory  does  not  account  for  shocks,  the 
disagreement  in  this  region  is  not  surprising. 
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The  next  wing  subjected  to  the  theoretical  analysis  was  a "warped" 
trapezoidal  wing  (See  figs.  3 and  10).  For  both  M » 1.61  and  2.01,  the 
leading  edge  Is  supersonic  and  multiple  "fold  over"  regions  occur  because 
of  the  low  aspect  ratio  and  the  fact  that  the  wing  Is  trapezoidally 
shaped.  The  equation  governing  the  mean  warped  camber  surface  Is 


(z/ c) * .033454  sin  (3irn/2)1  cos  (vOj  (38) 

Results  for  this  wing  are  presented  in  Figs.  15  and  16.  Considering 
the  wing  deformation,  the  agreement  is  remarkably  good  over  the  entire 
wing  for  both  Mach  numbers.  Note  also  that  the  addition  of  vorticlty 
paneling  to  the  ACp  distribution  predicts  the  right  trends  and  correctly 
accounts  for  local  wing  deformations  in  the  upstream  running  Mach  cones. 
Gapped  Elevon-Thickness 


For  supersonic  flow,  the  gapped  elevon  case  was  also  studied  in  the 


present  effort.  Little  experimental  data  is  available  for  comparison 
purposes;  however,  ore  previous  effort6  did  study  this  particular  problem 


using  a vortex  lattice  approach.  Conclusions  were  drawn  that  the  dis- 
agreement between  experiment  and  theory  was  attributed  to  viscous  effects. 
It  is  the  contention  of  the  authors  of  the  present  report  that  disagreement 


between  theory  and  experiment  is  due  not  only  to  viscosity,  but  also  to 
thickness  effects  and  blockage  of  part  of  the  elevon  immersed  in  the  wake 
of  the  "thick"  wing. 

To  validate  this  contention,  the  missile  configuration  used  in  Ref.  6 


was  analyzed  with  the  present  theoretical  approach  which  was  modified  for 
the  thick  wing-elevon  configuration.  A general  configuration  which  includes 


I 


a wing-body  and  full  span  elevon  has  not  as  yet  been  completely  analyzed 
and  hence  requires  additional  research. 

The  missile  configuration  under  study  consists  of  a slender  circular 
body  with  unwarped,  fixed,  low  aspect  ratio,  rear  mounted,  cropped  delta 
wings  and  constant  chord  elevons  (Fig.  3,  wing  #2).  Since  the  configura- 
tion operates  at  high  Reynolds  number,  viscosity  should  not  dominate  the 
pressure  field  so  that  potential  flow  theory  should  be  applicable.  Since 
both  the  wing  and  elevon  are  thin,  the  effects  of  thickness  on  the  lift 
are  negligible  and  both  may  be  treated  as  planar  lifting  surfaces. 

For  this  particular  analysis,  body  interference  effects  are  accounted 
for  with  slender  body  lift  ratios  as 


C„  * C [K  ,„.a  ] 

".<«)  \ W<B)  ” 

w 

where  C is  the  normal  force  on  the  wing  in  the  presence  of  the  body 

w(B) 

and  is  t*ie  lift  ratio  factor  for  the  wing-body  at  angle  of  attack. 

The  wing  is  assumed  to  lie  along  the  body  centerline.  Lift  ratios  are 
based  on  total  loads  so  only  the  net  effects  of  the  body  on  the  wing  or 


elevon  can  be  determined. 


The  elevon  is  located  immediately  aft  of  the  wing  as  shown  in  Fig.  17. 
The  flow  just  aft  of  the  wing  trailing  edge  is  turbulent,  as  with  any 
blunt  trailing  edge  and  impacts  on  the  elevon.  The  details  of  the  flow 
field  are  not  well  known  and  are  consequently  difficult  to  handle.  However, 


for  purposes  of  this  analysis,  it  was  assumed  that  the  wake  pressure  coef- 
ficient is  zero  such  that  essentially  freestream  conditions  exist  where 
the  wake  impacted  on  the  elevon.  In  short,  part  of  the  elevon  was 


r 


"blanked  out"  from  the  rest  of  the  flow  field  by  the  wing  wake.  Thus  the 
elevon  was  broken  up  into  two  separate,  smaller  lifting  surfaces,  and 
E 2 as  depicted  on  Fig.  17. 

The  flow  was  assumed  to  be  tangent  to  the  surface  of  the  aft  portion 
of  the  wing  for  all  angles  of  attack.  Thus  the  size  and  location  of  the 
"blanked  region"  on  the  elevon  was  the  same  for  all  the  no-gap  cases. 

When  the  wing  is  at  an  angle  of  attack,  the  downwash  it  creates 
reduces  the  angle  of  attack  of  the  flow  that  the  elevon  experiences.  No 
wing  wake  roll-up  occurs  for  this  case  and  thus  the  wake  of  the  wing 
downwash  was  assumed  to  be  flat.  Elliptic  wing  loading  and  angles  of 
attack  were  assumed  so  the  normal  force  and  lift  were  approximately  equal. 

The  downwash  angle  created  by  the  wing  in  the  presence  of  the  body  is 


given  as 


. = _ w(B) 

w it  AR 

w 


The  net  elevon  angle  of  attack  is  then 


a = a + e . (41) 

E w w 

The  angle  of  attack  contributions  of  the  elevon  characteristics  are  distinct 
from  the  deflection  angle  contributions.  However,  since  linear  aerodynamics 
is  assumed,  the  two  contributions  can  be  superimposed  to  determine  the 
overall  elevon  characteristics. 

The  strip  of  the  elevon  forward  of  the  "blanked  region,"  E^,  also 
generates  downwash  as  it  is  assumed  to  be  a separate  airfoil.  However, 
calculations  indicate  the  induced  angle  to  be  negligibly  small  in  com- 


parison to  the  elevon  deflection  angle. 

The  elevon  normal  force  and  hinge  moment  equations  with  the  thickness 


effects  included,  are 
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W V ‘ + +c»6  1*« 

v ' E E2 
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KE(B)  CH  aE  + kE(B) [CH{  + CH6  ]6E 
“E  Ex  E, 


(42) 


In  terms  of  the  wing  angle  of  attack  and  downwash  angle, 


W *««>  + ^ + V ,S* 

E E1  E2 


W WV'VV'WV'V*1 


(43) 


The  first  comparison  of  the  theory  with  experimental  data  is  shown 
in  Fig.  18.  For  this  case,  both  the  wing  and  elevon  are  modeled  as  zero 
thickness  flat  plates.  The  analytical  results  fall  below  the  experimental 
values  but  the  theory  does  give  the  correct  trends  for  normal  loading 
and  hinge  moment.  Figure  19  presents  the  results  for  the  wing  thickness 
effects.  The  angle  of  attack  trends  are  the  same  as  in  Fig.  18.  Now, 
however,  the  analytical  results  agree  with  the  experimental  results  much 
better,  particularly  the  normal  force  coefficients.  Wing  thickness, 
therefore,  exerts  substantial  influence  on  the  elevon  aerodynamic  charac- 
teristics, The  discrepancies  in  magnitude  and  slope  indicate  that  the  very 
simple  flow  model  used  here  is  not  completely  accurate. 

An  investigation  of  possible  wing  downwash  effects  on  the  size  and 


location  of  the  "blanked  region"  was  made.  For  this  model  it  was  assumed 
that  the  induced  downwash  at  the  wing  trailing  edge  caused  the  blanked  out 
region  to  shift  locations  on  the  elevon.  The  results  are  presented  in 
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Fig.  20.  The  poor  agreement  with  experimental  data  indicates  that  the 
assumed  flow  model  is  not  valid  and  that  physically  changes  in  angle  of 
attack,  i.e.,  changes  in  downwash,  do  not  cause  significant  alteration  in 
the  blanked  out  region  on  the  elevon. 

As  a means  of  reducing  the  elevon  hinge  moment,  a four-inch  cut  was 
made  into  the  trailing  edge  of  the  wing  so  that  a gap  was  formed  between 
the  wing  and  the  elevon.  In  this  instance,  the  elevon  was  treated  as  if 
it  were  not  "blanked  out."  A comparison  with  the  no-gap  configuration  is 
shown  in  Fig.  21.  The  theory  predicts  the  experimental  results  well,  if 
somewhat  conservatively.  Thus  the  wing  thickness  effects  are  only  sub- 
stantial if  the  elevon  is  very  close  to  the  wing  trailing  edge.  Addition- 
ally, it  was  observed  by  Selden6  that  the  elevon  behaves  as  if  it  were  in 
freestream  flow  when  it  is  deflected  -20°  or  more.  Therefore,  the  wing 
thickness  effects  appear  to  be  important  only  when  the  elevon  is  immediately 
aft  of  the  wing  and  when  it  is  slightly  deflected. 

In  conclusion,  it  is  quite  apparent  that  wing  trailing  edge  thickness 
does  indeed  have  significant  influence  on  the  loads  experienced  by  the 
elevon.  However,  it  should  be  noted  that  these  computations  are  only 
preliminary  and  that  a complete  model  of  the  flow  field  has  not  been 
established.  Further  research  using  numerical  techniques  (non-potential 
flow)  may  be  required  before  the  problem  is  finally  solved. 
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